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I. INTRODUCTION 



It has been more than a quarter century since Hawking's remarkable discovery |]J that 
a black hole is not completely black but can emit radiation from its event horizon. An 
important subject in black hole physics is to reveal the thermal properties of various black 
holes. Last few decades has witnessed much progress on investigating the thermal radiation 
of scalar fields or Dirac particles in some stationary axisymmetry black holes. Nevertheless 
most of these efforts (see f||-§J, for examples) were concentrated on studying the thermal 
properties of static or stationary black holes. Because a realistic black hole in astrophysics 
can radiate or absorb matter surrounding it, it is nonstationary and evolves in the time. 
Thus the thermal properties of non-stationary spacetimes are more interesting than that of 
static or stationary black holes, and worth much more studies about them. A well-known 
method to determine the location and the temperature of the event horizon of a dynamical 
black hole is to calculate vacuum expectation value of the renormalized energy momentum 
tensor |J. But this method is very complicated, it gives only an approximate value of the 
location and that of the temperature. Thus it is of limited use and meets great difficulties 
in many cases. 

To study the Hawking evaporation of the non-stationary black holes, Zhao and Dai || 
suggested a new method of the generalized tortoise coordinate transformation (GTCT) which 
can give simultaneously the exact values of the location and the temperature of the event 
horizon of a non-stationary black hole. By generalizing the common tortoise-type coordinate 
r* = r + ^- ln(r — r#) in a static or stationary spacetime 0,[7| (where n is the surface gravity 
of the studied event horizon) to a similar form in a non-static or non- stationary spacetime 
and allowing the location of the event horizon r# to be a function of the advanced time 
v = t+r* and/or the angles 8, (p, the GTCT method reduces Klein-Gordon or Dirac equation 
in a known black hole spacetime to a standard wave equation near the event horizon. For 
instances, the location of the event horizon is a constant (r# = 2M) in the Schwarzschild 
black hole while it is a function of the advanced time (r# = r#(t>)) in a Vaidya-type space- 
time. This method has been applied to investigate the thermal radiation of scalar particles 
in the non-uniformly accelerating Kinnersley "photon rocket" solution || and in the non- 
uniformly accelerating Kerr black hole |J as well. 
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However, it is very difficult to investigate the quantum thermal effect of Dirac particles in 
the non- stationary black holes. The difficulty lies in the non-separability of variables for the 
Chandrasekhar-Dirac equation fH| in the most general spacetimes. The Hawking radiation 
of Dirac particles has been studied only in some non-static black holes [fLl}| . Recently we 



T2| have tackled with the evaporation of Dirac particles in a non-stationary axisymmetric 



black hole. Making use of the GTCT method, we consider the asymptotic behaviors of 
the first-order and second-order forms of Dirac equation near the event horizon. Using the 
relations between the first-order derivatives of Dirac spinorial components, we eliminate the 
crossing-terms of the first-order derivatives in the second-order equations and recast each 
equation to a standard wave equation near the event horizon. Not only can we re-derive 
all results obtained by others ||13|| , but also we find that the Fermionic spectrum of Dirac 
particles displays another new effect dependent on the interaction between the spin of Dirac 
particles and the angular momentum of black holes. This spin-rotation effect is absent from 
the Bosonic spectrum of Klein-Gordon particles. 

It is natural to see whether or not our method can work effective in other cases. In this 
paper, we apply it to deal with the Hawking effect of Dirac particles in a non-spherically 
symmetric and non-stationary Kinnersley black hole, namely Kinnersley "photon rocket" 
solution ]TJ|]. The local event horizon of a dynamical black hole is determined here by the 
null hypersurface condition. We use different methods to deduce the location of the event 
horizon and find that they all give the same result. By means of a GTCT, we can also 
derive the event horizon equation from the limiting form of the first-order Dirac equation 
near the event horizon. The location and the shape of the Kinnersley black hole is not 
spherically symmetric ||. Then we turn to the second-order Dirac equation. With the 
aid of a GTCT, we adjust the temperature parameter in order that each component of 
Dirac spinors satisfies a simple wave equation after being taken limits approaching the event 
horizon. We demonstrate that both the shape and the temperature of the event horizon of 
Kinnersley black hole depend on not only the time, but also on the angles. The location and 
the temperature coincide with those obtained by investigating the Hawking effect of Klein- 
Gordon particles in the accelerating Kinnersley black hole JJ. But the thermal radiation 
spectrum of Dirac particles shows a new effect dependent on the interaction between the 



3 



spin of Dirac particles and the angular acceleration of black holes. This effect displayed 
in the Fermi-Dirac spectrum is absent in the Bose-Einstein distribution of Klein-Gordon 
particles. We find that this spin-acceleration coupling effect does not exist in a non-uniformly 
rectilinearly accelerating Kinnersley black hole. 

The paper is outlined as follows: In Sec. 2, we introduce the most general GTCT and 
derive the equation that determines the location of the event horizon from the null surface 
condition. Sec. 3 is devoted to discussing the Hawking radiation of Dirac particles in 
the Kinnersley spacetime. First, we work out the explicit form of Dirac equation in the 
Newman-Penrose |nj formalism and investigate the asymptotic behavior of the first-order 
Dirac equation near the event horizon. The equation that determines the location of the 
event horizon can be inferred from the vanishing determinant of the coefficients of the first- 
order derivative terms. Next, we use the relations between the first-order derivative terms to 
eliminate the crossing-term of the first-order derivatives in the second-order Dirac equation 
near the event horizon, and adjust the parameter k introduced in the GTCT so as to recast 
each second-order equation into a standard wave equation near the event horizon. In the 
meantime, we can get an exact expression of the Hawking temperature. Then the second- 
order equation is manipulated by separation of variables and the thermal radiation spectrum 
of Dirac particles are obtained by Damour-Ruffini-Sannan's method 0,0. In Sec. 4, we give 
a brief discussion about the new effect which represents the interaction between the spin of 
particles with spin-1/2 and the angular acceleration of black holes. 



II. KINNERSLEY BLACK HOLE AND ITS EVENT HORIZON 



The Kinnersley metric |fL4[|, generally called as the "photon rocket" solution, is inter- 



preted as the external gravitational field of an arbitrary accelerating mass. In the advanced 
Eddington-Finkelstein coordinate system [v;r;9;ip], the line element of Kinnersley's rocket 
solution reads 

ds 2 = 2dv(Gdv - dr - fr 2 d6 + gr 2 sin 2 Odkp) - r 2 (d9 2 + sin 2 9d<p 2 ) , (1) 

where 2G = 1 - 2M(v)/r - 2arcos9 - r 2 W*W, W = f - igsinO, W* = f + ig sin 9, 
f = b sirup + ccosip — a sin 9, and g = (bcosip — csimp) cot 9. The metric ([!]) is of type-D 



under Petrov classification. The arbitrary function M(v) describes the change in the mass 
of the source as a function of the advanced time; a = a(v), b = b(v) and c = c(v) are 
acceleration parameters: a is the magnitude of acceleration, b and c are the rates of change 
of its direction. The co-moving spherical coordinate system is oriented in such a way that 
the direction 9 = to the north pole always coincides with the direction of the acceleration. 

The spacetime geometry of an evaporating black hole is characterized by three surfaces: 
the timelike limit surface, the apparent horizon and the event horizon. According to York 



16], the horizons of a dynamical black hole may be obtained to first order in luminosity 
by note that (i) the apparent horizons are the outermost "trapped" surfaces that the ex- 
pansion of null-geodesic congruences (or null rays parametrized by v ) $ ~ 0; (ii) the event 
horizons are null surfaces where the acceleration of null-geodesic congruences d 2 r/dv 2 ~ 0, 
or equivalently they are determined via the Raychadhuri equation by the requirement that 
dd/dv m as they must be strictly null, and (iii) the timelike limit surfaces are defined as 
surfaces such that g vv = 0. 

It is generally accepted that the event horizon is necessarily a null surface and is defined 
by the outermost locus traced by outgoing photons that can "never" reach arbitrarily large 
distances fl(|. In a nonstationary black hole spacetime, the event horizon should still be a 



null surface that satisfies the null surface condition: g^diFdjF = 0, and the event horizon 
determined by the above null hypersurface condition is, in fact, a local event horizon. We 
shall adopt this definition and use different methods to derive the equation that determines 
the location of local event horizon of an arbitrarily accelerating Kinnersley black hole. We 
find that each method can give the same result consistently. 

First, let's seek the local event horizon of Kinnersley spacetime (|1|) by using of the null 
surface condition. From the null surface equation F(v,r,6,ip) = 0, namely r = r(v,8,(p), 
one can easily obtain 

d v F + d r Fd v r = , d e F + d r Fd e r = , d v F + d r Fd^r = . (2) 

Substituting into the explicit expression of the null surface condition g^diFdjF = in 
the Kinnersley metric (|l|) 

(2G + r 2 W*W)(d r F) 2 + 2d r F(d v F - fd e F + gd v F) 



1 - -x 2 1 



+ -{d e FY + -^—^—(d^F) = , (3) 
one gets 



2G + r 2 W*W - 2d v r + 2fd e r - 2gd lf r + + ff^ n = . 

The local event horizon is the hypersurface r = rn{v, 8, <p) that satisfies the above equation 
or 

2M 

1 2ar H cos 9 - 2r H v + 2fr H e - 2gr H ip 

r H 

in which r HjV = d v r H , r H)d = der H and r H>lfi = d^ru can be viewed as parameters depicting 
the evolution of the event horizon. 

When a = b = c = but M/fl, the event horizon of Vaidya black hole is located at 
rn = 2M/(1 — 2rfj, v )] When M = b = c = and a = const, the Rindler event horizon of 
the uniformly rectilinearly accelerating observer satisfies 

1 — 2arn cos 6 — 2aru e sin 6 H J- = , 

it is a paraboloid of revolution m = l/a(cos# ± 1). 

In the case where M(v), a(v), b{y) and c(v) are not equal to zero, Eq. (|) has in 
general three roots. In this general case, the analysis is a little involved, and will not be 
discussed here. There should exist two kinds of event horizon: Rindler-type horizon and 
Schwarzschild-type horizon. All of them depend not only on v, but also on 9, (p. It means 
that the location of the event horizon and the shape of the black hole change with time. 
The location of event horizon is in accord with that obtained in the case of discussing about 
the thermal effect of Klein-Gordon particles in the same spacetime ||. 

Next, we adopt the GTCT method to deduce the equation of local event horizon. As 
the Kinnersley metric is lack of any symmetry, we introduce the most general form of the 
GTCT § as follows 

r* = r + - 1 - In (r - r H ) , v m = v — v , 
2k 

6* = 8 - 9 , (f* = (p - (fo , (5) 
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namely, 



dr* = dr + — (dr - r Hv dv - r H0 d9 - r H ^dip] , 

2k[t — r H ) v ' 

dv* = dv , d9* = d9 , d(fi* = dip , 

where r^ = rn{v,9,ip) is the location of the event horizon, k is an adjustable parameter 
and is unchanged under tortoise transformation. All parameters vo, 9o and <po are arbitrary- 
constants characterizing the initial state of the hole. 

Applying the GTCT ([5]) to the null hypersurface equation (||) and taking the r — > 
r H(v ,9 ,{p ), v — > v , 6 — > and — > <po limits, the event horizon equation is then 
obtained by letting the term in the bracket before (^-F) 2 to be zero, [] 

2 2 

2G - 2r H , v + r 2 H W*W + 2fr Hfi - 2gr H ^ + ^ + J H % D = . (6) 

r A H r H sin Uq 

Eq. (||) is just the same equation (|]) when v = v , 9 = 9 and y? = ip . Because here we deal 
with the case of a slow evaporation of black holes, we need later only consider the situation 
very close to the initial state of the event horizon, namely r# w i"h(vo, 9q, (p ) when v ~ Vq, 
9 w 9 and ~ <^o- This assertion is due to that the GTCT approach is a local analysis 
method, the latter originating from that Hawking radiation comes from vacuum fluctuations 
near the event horizon. 



III. HAWKING EVAPORATION OF DIRAC PARTICLES 

Now we turn to investigating the quantum feature of the Kinnersley spacetime, especially 
the thermal radiation of electrons, that is, we must derive the Hawking temperature of the 
event horizon and the thermal radiation spectrum of Dirac particles from the event horizon. 
To this end, we work out the spinor form of Dirac equation in the Newman- Penrose (NP) 



Throughout the paper, we make a convention that all coefficients in the front of each derivatives 
term take values at the event horizon ru = rn(vo, @o, y?o) when a GTCT is made and followed by 
taking limits approaching the event horizon, e.g., G, /, g and W take their values at v = vq, 9 = 9q 
and ip = ifQ. 

7 



15 1 formalism. We choose a complex orthogonal null-tetrad system in the Kinnersley black 
hole such that its directional derivatives are 



D 



-d r 



A = a. + Gd r 



<) = - r 2 Wd r + d e + -^—d, 

■ \ sin H 



V2r 
1 



5 = -r 2 W*d r + d e 



sinf 7 
sin 



(7) 



It is not difficult to determine the non-vanishing complex NP spin coefficients in the above 
null-tetrad as follows f\ 





1 


G 


p = 




fi = — 




r 


r 




W 




T = 




7T = - 


V = 


1 

V2r 


'{2rG - 



-G r + ig cos6 l )/2 



W* 
72 



a 



cot e w* 
+ 



2V2r V2 



P 



cote 



(2rG - r 2 G, r )W* + r 2 W* v + G, e 



iG, 



sin 6 1 



The dynamical behavior of spin- 1/2 particles in curved spacetime is described by the 
four coupled Chandrasekhar-Dirac equations [[HJ expressed in the following spinor form 

(-D + e — p)F\ + (5 + tt — a)F 2 = % ^G 1 , 

V2 

(A + / x- 7 )F 2 + (5 + /3-r)F 1 = ^G 2 , 

(£> + e* - p*)G 2 -(5 + r- a*)G l = ^F 2 , 

v2 



(A + p* - -f)G 1 - {6 + (3* - r*)G 2 



V2 



Ft 



(9) 



where po is the mass of Dirac particles. Inserting the needed NP spin coefficients and making 
substitutions Pi = \/2rFi, P 2 = F 2 , Qi = Gi, Q2 = \[2rG 2 into Eq. @, we obtain 

-V P 1 + (£ - r 2 W*V 2 )P 2 = z f i rQ 1 , 
2r 2 B\P 2 + (£ f - r 2 WV )P 1 = ip rQ 2 , 
-V Q 2 - (£t _ ^WX^Qi = W P 2 , 

2r 2 B 1 Q 1 - (£ - r 2 W*D )<2 2 = W^i , (10) 



2 Here and hereafter, we denote G r = dG/dr, etc. 
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in which we have defined operators 



Bn = d v + GV n + (G r - ig cos 9) /2 



£ = d e + - cot 9 
2 



i 



d r + n/r . 
i 



^—;d v , £ f = d e + ^ cot 9 + — 7A 
sinfl 2 sm9 



Because the Chandrasekhar-Dirac equation (^) can be satisfied by identifying Gi, G 2 with 
F 2 *, —F*, respectively, so one may deal with a pair of components (Pi,P2) only. Although 
Eq. ([K]) can not be decoupled, to deal with the Hawking radiation, one should be concerned 
about the asymptotic behavior of Eq. ( |10D near the horizon only. Under the transformation 
Eq. (|II]) with respect to (Pi, P2) can be reduced to the following limiting form near the 
event horizon 







-z-Pi + (r H ,6 
or. 



I r _|_ r 2 _^_p 2 

sin 9 ' 9 H / dr* 







(rn,e + -^-J^ + rlW) ±P, - 2r\ (G - r H>v ) ^-P 2 







11) 



dr* 

after being taken the r — > r H (v ,9 ,(p ), v — > v , 9 — > 9 and ip — > (p limits. If the 
derivatives -§~P\ and -§~P% m Eq. (|ll]) are not equal to zero, the existence condition of 
non-trivial solutions for Pi and P2 is that its determinant vanishes, which gives the above- 
head equation (^|). This treatment can be thought of as another derivation of the location 
of event horizon. It is interesting to note that a similar procedure applying to another pair 
of components [Qx^Q^) will bring about the same result. 

To investigate the Hawking radiation of spin-1/2 particles, we need deal with the behavior 
of the second-order Dirac equations near the event horizon. A direct calculation gives the 
second-order form of Dirac equations for (Pi, P2) components as follows 



2r 2 B V + (£ - t^W^V-^ip - r 2 WV ) - /2 2 r 2 ]P 1 



+2r 2 (£ - r 2 W*V x )B\ - B (C - r 2 W*V 2 ) 



P. = 



(12) 



and 



2r 2 £>i#I + (£ f - ^WV^iC - r 2 W*V 2 ) - /J, 2 r 2 \P 2 
+ [£>_i(£ t ~r 2 WV ) - (£ f -r 2 WV x )V^P x =0. 



(13) 
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and 



Given the GTCT in Eq. ([5]) and after some lengthy calculations, the limiting form of 
Eqs. ( P^pTSD , when r approaches r H (v , 9 , ip ), v goes to v , 6 goes to 9 and tp goes to <p , 
yields 



/CPi 



A + r 2 H G tr + 2r 2 H G + 2r%W*W + 2ig cos 6 - 2r 2 H f cot 6 



+ cot 8 

i 



TH,e + 



sin Oc 



T H fi8 + 



vV 



sin O <9r* 







Pi + 2r 2 H 



G, 



sin #o 



G> + r 2 ^; + W\2r H G - r 2 H G, r ) 



_d_ 

dr* 



Po = 



(14) 



JCPo 



A + 3r 2 H G :r + r H (6G - Ar H , v ) + Qr%W*W - 2r 2 H f cot 9 



Agr H + i cot 9 ) - 



sin 9 



rH. 



(rH,ee + . ■■> 
v sin 6/, 



<9r* 



P 2 = 0. 



(15) 



In the above, the operator stands for the term involving the second-order derivatives 



A 
2k 



2r 2 H (2G - r HjV ) + 2r H W*W + 2r 2 H {fr Hfi - gr H; , 



2 



2-(fr H + r Hi( 



d 2 



+ 2( 



or^dv* v ' Or^dQ* v sm z 6*0 7 or*dp* 



d 2 



By the event horizon equation (§) or (|6|), we know that the coefficient A is an infinite limit 
of 0/0-type. Using the L' Hospital rule, we arrive at its obvious result 



A = lim 

r—*rn 



2r 2 (G - r H ,v) + r 4 W*W + 2r\fr Hfi - gr H)V ) + r% e + ^ 



r — Th 

2r%G >r + 2r%W*W - 2r u \r 2 Hfi + r 2 H J sin 2 9 ) . 



Now let us select the adjustable parameter k in the operator /C such that 
r H = ^~ + 1r\(2G - r H>v ) + 2r H W*W + 2r 2 H (fr H , e - gr H>v ) , 



(16) 



which means the surface gravity of the horizon is 



K 



r 2 H G, r + r H W*W - r H \r 2 Hfi + r%J sin 2 fl ) 



(17) 



r|(l - 2G) - r H W*W + r\ fi + r\J sin 2 6 ' 

With such a parameter adjustment and using relations (|TT|), Eqs. (|l^ , |i"5|) can be recast into 
the following standard wave equation near the horizon in an united form 
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-d 2 2 d 2 2C d 2 2Q d 2 d I 

.^ r 2 dr*dv* 3 ' dr*d9* dr*dip* 2 1 <9rv 

where f2, C3, C2, Ci shall be regarded as finite real constants. The physical meaning of Q is 

that it can be interpreted as the angular velocity of the black hole due to deformation. For 

completeness, they are listed as follows 



th, v „ , , r H . 

r 2 H sin 2 9 r H 



V = g- 3 In > c 3 = f + 



2 ' 



while both C\ and C2 are real, 



2G cot 6 / ir g , y 

G, r + 2igcos9 3— (r ff) e r-f- 

r# r H sin t^o 



+ 



C-? — sin ft 



2/cot0 o + — [rke + -^r) ~ -{rH.ee + -^r) (19) 



(G - r H)V )r 2 H 



sin 9 ' r H ^ ' sin 9 
r 2 H W* v + W*(2r H G - r 2 H G, r ) + G >6 - -^-G )V 

Sill t7n 



for = Pi, and 



2C 

2(C 2 + = — +G r + 2r H W*W - 2f cot 9 
r H 

cot / ir H , v s 1 / rff iW \ 

2~ r ^ + ^^r T\ r H,ee+ . 2n ( 20 ) 

r# v sm^ y r 2 H \ sin 9 Q J 

for if? = P 2 . The expression in Eq. ( [I9| ) is very complicated, but one can notice that the last 

term in the square bracket is proportional to u(th) (i.e. the value of the spin coefficient v 

at the event horizon r#). 

Because all coefficients in Eq. (|IBD can be viewed as constants approximately, the wave 

equation can be manipulated like an ordinary differential one. Now separating variables as 

\fj = R( r *)<d(9*)$(ip*)e- iuJV * and substituting it into Eq. ©, one gets 



0' = A6, $' = ((7 + im)$, 

R" = 2i(oo -mn-Ci + iC )R' , (21) 

where Co = C2 — AC3 + aQ, A is a real constant introduced in the separation variables, u 
the energy of electrons, m the quantum number of its azimuthal angular momentum. The 
solutions are 

R = flie 2*(w-mn-Ci)r.-2C r. + Rq _ ( 2 2) 
11 



The ingoing wave and the outgoing wave to Eq. fll8|) are, respectively, 

* in — B , 

*out = ^ in e 2i ^~ mn - Cl ) r *- 2Cor * , (r > r H ) . (23) 

Near the event horizon, we have r* ~ (2k) _1 ln(r — ?"•#)■ Clearly, the outgoing wave ^ ut(^ > 
r#) has a logarithm singular and is not analytic at the event horizon r = r#, but can be 
analytically extended from the outside of the hole into the inside of the hole through the 
lower complex r-plane (i.e., (r — r#) — > (rj? — r)e~ t7T ) to 

According to the method suggested by Damour and Ruffini |2j and developed by Sannan 
0, the relative scattering probability of the outgoing wave at the horizon is easily obtained 

^cmt 2 __ e -27r(w-mft-C*i)/K (25) 

The thermal radiation spectrum of Dirac particles from the event horizon of the hole is given 
by the Fermionic distribution 

M")> = e^-JtL + x » ( 26 ) 

in which T(u) is the barrier factor of certain modes, and the Hawking temperature Tjj = 

k/(2tt) is obviously expressed as 

T = 1 x Mr H - r 3 H a cos 9 - r 2 Hfi - r 2 H J sin 2 9 Q 
H Aitr H X Mr H + r 3 H a cos 9 + 2^{r 2 H>6 + r 2 H J sin 2 9 ) ' 1 J 

It follows that the temperature depends not only on the time v, but also on the angles 9 

and ip because it is determined by the surface gravity k, a function of v, 9, and ip. The 

temperature is consistent with that derived from the investigation of the thermal radiation 

of Klein-Gordon particles in the Kinnersley black hole M. 



IV. SPIN-ACCELERATION COUPLING EFFECT 

There are two parts in the thermal radiation spectrum (|26|), one is the rotational energy 
mfl; another is C\ due to the coupling between the spin of electrons and the angular mo- 
mentum of the black hole. Comparing the thermal spectrum of spin-1/2 particles with that 

12 



of scalar particles ||, we find that an extra term C\ appears in the former and is absent in 
the latter. From its explicit expression, 



C 3 sin 6n 



G, 



C 1 = \ r 2 H g, v + (2r H G - r 2 H G, r )g - —g- 

sin On 



2(G - r HA 
fl sin 9n 
'2(G-r HtV ) 



rj I f jV + {2r H G-rj I G, r )f + G jl 



+°™°° + %0T o - <* = «> <* 

r Hftp cos go 
2r 2 H sin 2 9 



Ci = -S^S7 > (* = ^) (29) 



one can find that G\ — for \I/ = Pi, P2 in the non- uniformly rectilinearly accelerating black 
hole (b = c = and = 0). Further more, if we only consider the second term in Eq. 
(|28|), then we can rewrite the "spin-dependent" term as 

uj s ~ gr f' y T/° , (a = 1/2, * = Pl ; s = -1/2, v& = P 2 ) (30) 
sin u 

for different spin states s = ±1/2 respectively. Here we are interested especially in this term 
because it is obviously related to the spin of electrons in different helicity states. The factor 
r H,ip describes the deformation of black hole during its evolution, while the factor cos 9q comes 
from the scalar product between the spin vector of electrons and the "deformation angular 
momentum" rjj tV / {r 2 H sin 2 9q) of the black hole. One can notice that this "deformation 
angular momentum" is also contained in the expression of Q. Thus this new term represents 
the spin-acceleration coupling effect which can be interpreted as the interaction between the 
angular variation of the black hole and the spin of the particles. 

One can observe that u s vanishes at the equator 6 = it/ 2 and diverges in s9~ 2 near the 
north pole 9 w and in — s(n — 9)~ 2 near the south pole 9 w ir. In general, the dependence 
of C\ upon the mass M and the acceleration parameters (a, b, c) is indirect and non-trivial 
because they are involved in the expression of r#. To understand this point better, we 
consider a particular case where we take M, a, b and c as constants and set 9 = 9 = ir/2. 
In this case, we have 2G = 1 — 2M/r — r 2 f 2 , f = bsimp + ccosy? — a and g = as well as 
r H,v = r H,e = 0. Thus the event horizon equation becomes 

r 2 H -2Mr H + r 2 H ^ = 0, (31) 
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and it has solutions r# = M(l ± cosy); M(l ± siny), which demonstrate that the event 
horizon is a triangle function of y. Without loss of generality, let's take one solution as 
th = M(l + cosy). Then the Hawking temperature is 

rp = cosy 

27rM(l + cos y) (3 — cos y) 

Further, we have C3 = and Q = siny/M(l + cosy) 2 . For \P = P 2 , G\ vanishes; while for 
^ — P\ due to f jV = G y e = 0, it becomes 

q sin y (cosy — 2) (6 sin y + ccosy — a) 

sin 2 y + M 2 (l + cosy) 4 (6siny + ccosy — a) 2 

Obviously G\ is linearly dependent on a, b and c if M = 0. We note that r# = 2M, 
T H = 1/(8ttM) and tt = C x = at y = y = 0. 

In the general case where M, a, b and c are arbitrary functions of the advance time v, 
to analysis Ci, Q and Tjj is, however, very complicated. It needs us first to determine the 
explicit expression of r# in terms of these parameters, but this topic apparently exceeds the 
context of this paper. We wish to discuss them in other circumstances. 



V. CONCLUSIONS 

Equations (||) and ( p7]) give the location and the temperature of the event horizon of 
the Kinnersley black hole, which depend not only on the advanced time v but also on 
the angles 9,ip. Eq. (|26|) shows the thermal radiation spectrum of Dirac particles in the 
arbitrarily accelerating Kinnersley spacetime, in which an extra term C\ appears. We find 
that C\ vanishes when b = c = and th, v = 0. This means that this term does not exist 
in the rectilinearly accelerating Kinnersley black hole. We contend that this new effect 
probably arise from the interaction between the spin of Dirac particles and the acceleration 
of the evaporating black hole. Besides, we notice that this spin-acceleration coupling effect 
appearing in the Fermionic thermal spectrum of Dirac particles is not shown in the Bosonic 
spectrum of scalar particles. 

In summary, we have studied the Hawking radiation of Dirac particles in an arbitrarily 
accelerating Kinnersley black hole whose mass changes with time. The location and the 
temperature of the event horizon of the accelerating Kinnersley black hole are just the same 
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as those obtained in the discussion on thermal radiation of Klein-Gordon particles in the 
same spacetime. But the thermal spectrum of Fermi-Dirac distribution of particles with 
spin-1/2 displays an extra interaction effect between the spin of Dirac particles and the 
angular acceleration of black holes. The character of this term is its obvious dependence on 
the different spin states. This spin-acceleration coupling effect is absent from the thermal 
radiation spectrum of Klein-Gordon particles. In addition, our discussion here is easily 
generalized to the most general case of a non-uniformly accelerating Kinnersley black hole 
with electric charge Q(v), magnetic charge P(v) and cosmological constant A (2G — 1 — 
2M(v)/r + (Q 2 + P 2 )/r 2 - 4a cos 6(Q 2 + P 2 )/r - 2ar cos 9 - r 2 (f 2 + g 2 sin 2 9) - Ar 4 /3). 
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